A dimerized spin fluid in a one-dimensional electron system 
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Abstract 

The ground state of a one-dimensional Hubbard model with a bond-charge 

attraction W term at half-filling is investigated by the density matrix renor- 

malization group method. It is confirmed that the spin gap will be closed at 

U > 8W. But the long-range bond order wave survives even when the spin 

gap is closed. It indicates that the ground state is a novel dimerized spin fluid 

at U > 8W. By a charge-spin transformation, it is shown that there should 

be a dimerized metallic phase at U < —8W. Furthermore, it is found that the 

Hubbard interaction U enhances initially the dimerization for a weak bond 

charge attraction W whereas it reduces monotonously the dimerization for a 

stronger bond charge attraction W. 
PACS:71.30.+h; 71.10.Fd; 71.10.Pm 
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The response of correlated electron systems to lattice distortions has been extensively 
studied motivated not only by theoretical interest but also by the discovery of quasi- 
one-dimensional conductors and high-Tc superconductivity. It is well known that Peierls 
instabilityB could induce a lattice dimerization in a one-dimensional(ld) electro n-phonon (e- 
p) interacting system! at half filling. A dimerized system has a long-range bond order wave 
(BOW). The BOW could be enhanced by an electron-electron (e-e) interaction, although a 
strong Coulomb e-e interaction will always reduce the BOW, which exists until U ooi@. 
In the limit of very strong on-site Coulomb interaction, the Peierls instability is usually re- 
ferred to as the spin-Peierls instability. Within a spin-Peierls model, the quantum phonons 
have been taken into account recentlylHi and a quantum phase transition from a gapless 
spin-fluid state to a gapped dimerized phase is found to occur at a nonzero value of the 
spin (electron) -phonon coupling. It is believed that the BOW long-range order will disap- 
pear as the spin gap is closed!. On the other hand, the BOW has been shown to exist in an 
interacting electron model recently.!^ In this paper, we investigate the Id Hubbard model 
with a bond-charge attraction term at half-filling by the density matrix renormalization 
group (DMRG) method0 and find an interesting phenomenon that the BOW long-range 
order survives even the spin gap is closed for U > 8W, which shows that the ground state 
is a novel dimerized spin fluid phase. The result is in contradiction to that obtained by a 
weak-coupling continuum-limit approachlll. 

The Hamiltonian of the Id interacting electron system at half-filling we consider is 
as follows 

H = Hu + Hw , (1) 
where Hu is the usual Hubbard model 

Hu = -tYl {4,aCi+i,. + c;+i,.Q,,) - I) Ul - I) , (2) 



and Hw is a bond-charge attraction {W > 0) term 



Hw = -WY.iBi,i+,y , (3) 



with the bond-charge density operator defined as 

Bl,l + 1 = J2 (4,aC/ + l,<x + 4^1,aCl,a) ■ (4) 
a 

The bond-charge attraction Hw can be obtained from a Su-Schrieffer-Heeger e-p 
interaction^ by integrating out the phonon degree of freedom in the antiadiabatic hmit 
(M So the W-term could also be viewed as a contribution of the e-p interac- 
tion with the parameter related to the dimensionless e-p coupling^ by ly = TiXt/A and 
the Hamiltonian (1) has included both the quantum phonon fiuctuation and the electron 
correlation for a one-dimensional electron-phonon interacting system. 

First, let us consider the symmetry properties of the model (1). It is obvious that the 
model has a usual SU(2) rotational symmetry in spin space. The generators of the rotations 



are 



i i i 

which commute with the Hamiltonian (1). There is also an SU(2) axial charge symmetry as 
in the Kondo mo delS and the Hubbard mo deM The axial charge generators are 

= ^ E (4tc.t + iic^,i - 1) , = E (-1)^ ^T^a, = E i-^y (6) 

i i i 

which have the same commutation relation as the spin operators in Eq.(5). So the Hamil- 
tonian (1) has an SU(2)xSU(2) symmetry group as the alternating Hubbard modeS It is 
clear that the z component of the axial charge operator is nothing but the total electron 
number operator, while its transverse components are the staggered pairing operators, called 
?7-pairing.ii 

Due to the commutation between the spin and charge generators, we could identify the 
eigenstates of the Hamiltonian (1) by the eigenvalues of the total spin 5*^ and total charge 
operators P \ S^, P), then the spin gap and charge gap are defined as 

A, = Eo (1, 0) - Eo (0, 0) , A, = Eo (0, 1) - Eq (0, 0) (7) 
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where Eq {S^, P) is the lowest energy in the (S*^, P) subspace, and it must have the eigen- 
values for the spin S = and charge / = J^H. 
The charge-spin transformation 

Ci,T ^ Ci,|,Q,| ^ (8) 

interchanges the spin generator (5) and charge generator (6), hence the spin and charge gaps 
in Eq.(7). It is straightforward to show that the bond-charge operator (4) is invariant under 
the transformation that Hubbard repulsion turns into attraction . Thus we have an equality 
for the spin and charge gaps 

A,(t,W^,f/) = A,(t,W^,-t/) (9) 

So if the spin gap is closed at a critical [/ (> 0), the charge gap must be closed at — t/, the 
system will be in a metallic phase. 

Now, we come to the numerical treatment of the Hamiltonian (1) by the powerful DMRG 
method^, which allows an essentially exact treatment of electron correlations at a fully 
quantum mechanical level. The open boundary condition is adopted since it is usually 
performed much better than the periodic boundary condition. The hopping integral t is 
set to 1 as the unit of energies. Up to m = 512 states are tested and finally m = 128 
states are kept in the DMRG calculation for the systems with L = 48, 64, 80, and 128, the 
truncation error is less than 10~^ — 10~® at all steps. After four iterations, the change in 
the ground-state energy is at seven digits between the successive iterations. We checked our 
DMRG calculations against exact numerical results for long (up to 128 site) non- interacting 
{ U = W = 0) chains, excellent agreements were found both in the correlation functions 
and in the energies of the ground and low-lying excited states. The corresponding relation 
(9) between the spin and charge gaps is also confirmed in our calculation within numerical 
error. Moreover, we have carefully compared our charge gap with the Bethe Ansatz exact 
solution of the Hubbard model@ and obtained quantitative agreement. 

The first quantity we calculate is the spin gap, Eq. (7). In the case of U = 0, the bond- 
charge attraction W-term will induce a BOW long-range order as the e-p interaction in the 



SSH mo deE. The charge and spin gaps are equal in its excitation spectra. All fluctuations 
besides the BOW are completely suppressed. In the presence of Hubbard repulsion U > 0, 
the charge gap increases while the spin gap decreases with increasing U. Figure 1 shows the 
spin gap As versus the Hubbard U for different values of the bond-charge attraction W. The 
inset shows a linear extrapolate of Ag (L) with the inverse of chain length 1/ (L + 1). The 
spin gap As at thermodynamic limit in Fig. 1 decreases monotonically with the increase of 
U. 

Within a g-ology continuum-limit model, Emeryil has given out the phase diagram for 
an interacting fermion system in the weak coupling region by using the renormalization 
group method. The interaction parameters in the Hamiltonian (1) can be transformed into 
the g-ology ones, 

g^ = U-8W, g2 = U, gs = U + 8W. (10) 

The ground state depends only on the sign of gi. 5^1 = is a phase transition point from a 
charge density wave (CDW) to a spin density wave (SDW), where the BOW character was 
not considered. Combining Emery's theory and bosonization method, Japaridze considered 
the BOW order and confirmed the phase transition from a spin-gapped state to a gapless 
spin fiuid at f/ = 8W. Our calculation on the spin gap is in good agreement with that 
of the weak-coupling theoryHEl. The exact transition point is not easy to be identified 
numerically since the spin gap decreases exponentially, which is the same as that in the 
spin-Peierls modehlli 

We now investigate the nature of the BOW ordering by examining the staggered bond 
charge correlation function 

CbOW (r) = — {-ly^ {Bi^i+iBi+r,i+r+l) (11) 

^ i 

The spin density wave (SDW) order can be described by a similar staggered spin-spin cor- 
relation function 

CsDw{r) = \{-irY.{^>Ur) (12) 



where the z-component spin is defined as 

Si=^KT-^u)- (13) 

As is well known, for a uniformed system, the bond charge (Sj.j+i) is uniformly distributed 
along the chain, and a rapid decrease in the bond charge correlation function (11) is ex- 
pected, that is, the BOW order is short-ranged. However, for a dimerized system, the bond 
charge j+i) shows an oscillation on bonds, the staggered bond charge correlation func- 
tion Cbow i''^) "will approach to a nonzero constant AboW) which shows the existence of a 
BOW long-range order. Usually, it should take a much longer chain for the bond charge 
calculation to identify the BOW long-range orderlll, however, the bond charge correlation 
will be saturated for a chain with a hundred sites. 

In Figure 2 we show both the staggered bond charge (a) and spin-spin (b) correlation 
functions. From Fig. 2(a) we can see that the staggered BOW correlation functions show 
a similar behavior for both the spin-gapped state {U < 'SW ) and the gapless spin fluid 
state {U > 8W ). And at the same time, from Fig. 2(b), we can see that the spin-spin 
correlations show a different behavior in these two regimes, it has a short-range SDW order 
in the spin-gapped state ( U < 8W) and a quasi-long-range SDW order in the gapless spin 
fluid state {U > 8W). The line in Fig. 2(b) is a 1/r fltting. This result is in conflict with 
what Japaridze obtained in Ref. 15, where it was pointed out that the phase transition 
at [/ = 8W is from a spin-gapped state with a true long-range BOW order to a gapless 
spin-fluid state characterized by an identical 1/r decay for both quasi- long-range SDW and 
BOW orders. But our calculation indicates that the BOW correlation shows a long-range 
order even when the spin gap is closed in Fig. 2(a). While the obtained spin-spin correlation 
function can be fltted very well by a 1/r-decay function (the line in Fig. 2(b)), the bond 
charge correlation function in Fig. 2(a) obviously does not show the 1/r-decay behavior. It 
indicates that there is a novel dimerized spin fluid phase for the Hamiltonian (1) at U > 8W. 
By the charge-spin transformation (8), the ground state will be in a dimerized metallic phase 
for —U > 8W > 0, where the staggered BOW correlation shows a long-range order with a 
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vanished charge gap. 

Figure 3 shows the BOW order parameter Abow versus the Hubbard repulsion U for 
different bond charge attraction W for a L = 128 chain. In the weak bond charge interaction 
regime, the BOW order is first enhanced by the Hubbard repulsion U and reach a maximum, 
then the BOW order is reduced, such behavior is similar to the e-p model at adiabatic 
limit.ilS However, for a stronger bond charge interaction (such as W = 0.2), the BOW will 
be reduced monotonously by the Coulomb interaction. When the Hubbard U is increased 
further, we will see that the BOW order becomes almost independent of the bond charge 
interaction since the Hubbard U dominates. There another phase transition is expected for 
the BOW long-range order since the spin-Peierls model shows the transition!^ when the 
quantum fluctuations are taken into account. 

In summary, we have found a novel dimerized spin fluid phase for the Id Hubbard 
model with a bond charge attraction W-term. The result indicates that there are two phase 
transition points, the first one is around U = 8W for the spin gap and the second one is at 
t/ < oo for the BOW long-range order. The BOW long-range order is shown to exist even 
when the spin gap is closed when U > 8W. The charge-spin transformation gives that the 
ground state should be in a dimerized metallic phase for U > —8W. 

We would like to acknowledge the useful discussions with Prof. X. Sun, Y.P. Wang, Z.Y. 
Weng, T. Xiang, and G.M. Zhang. This work is supported by National Natural Science 
Foundation of China (No. 19725414) and partially supported by a grant from the Research 
Grants Council of the HKSAR, China (Project No. CUHK 4288/OOP). 
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FIGURES 

FIG. 1. The spin gap versus the Hubbard repulsion U for different values of the bond-charge 

attraction W. The inset shows a linear extrapolation of Ag (L) with the inverse of the chain length 
1/(L + 1). 

FIG. 2. (a) The staggered bond charge correlation function CBOwir), and (b) the staggered 
spin-spin correlation function CsDw{r) for the chain of length L = 128. The dashed line in (b) is 
al/r fit. W/t = 0.1. 

FIG. 3. The bond charge order parameter Abow versus the Hubbard repulsion U for the chain 
of length L = 128. 
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